Applicability of layered sine- Gordon models to layered superconductors: 
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In this paper, we propose a quantum field theoretical renormalization group approach to the vortex dynamics 
of magnetically coupled layered superconductors, to supplement our earlier investigations on the Josephson- 
coupled case. We construct a two-dimensional multi-layer sine-Gordon type model which we map onto a gas of 
topological excitations. With a special choice of the mass matrix for our field theoretical model, vortex domi- 
nated properties of magnetically coupled layered superconductors can be described. The well known interaction 
potentials of fractional flux vortices are consistently obtained from our field-theoretical analysis, and the physi- 
cal parameters (vortex fugacity and temperature parameter) are also identified. We analyse the phase structure of 
the multi-layer sine-Gordon model by a differential renormalization group method for the magnetically coupled 
case from first principles. The dependence of the transition temperature on the number of layers is found to be 
in agreement with known results based on other methods. 

PACS numbers: 74.20.-z, 74.25.Dw, ll.lO.Hi, ll.lO.Gh, ll.lOKk 

I. INTRODUCTION 

Recently, we have shown that layered sine-Gordon type models are probably not suitable for the description of Josephson- 
coupled layered superconductors, because the linear, confining potential that binds the vortices together cannot be obtained 
from the interaction of the topological excitations of the model, no matter how the interlayer interaction term is chosen iH]. 
On the other hand, vortex dominated properties of high layered superconductors and other types of layered materials^e.g. 
superconducting sandwiches, have already received a considerable amount of attention (see, e.g., Refs. 

[ISIHlII^IIlllSllllIMl), and the intuitively obvious connection of sine-Gordon models to these materials makes one wonder 
if at least, a restricted applicability of the layered, field-theoretical model persists. We also observe that recently, there is an 
increasing interest in the literature iflTl [Tsl [l9ll in constructing sine-Gordon type field theoretical models in order to understand 
better the vortex dynamics in layered superconducting systems. Our aim in this paper to follow this route by constructing a 
two-dimensional multi-layer sine-Gordon type model which can be used to describe the vortex behaviour of magnetically as 
opposed to Josephson-coupled layered superconductors, and to contrast and enhance our recent investigation |]J]. 

In a two-dimensional (2D) isolated superconducting thin film, the Pearl vortices are naturally identified as the topolog- 

ical excitations and can be considered as the c harg ed analogues of the vortices in a 2D superfluid which generate the Kosterlitz- 
Thouless-Berezinski (KTB) phase transition luOll . The logarithmic interaction between the vortices of the superfluid extends 
to infinity and as a consequence they remain bound below the finite KTB transition temperature (Iktb) dissociate above 
it ll20ll . Since the Pearl vortices carry electric charge, they always remain unbound due to the screening length Aoff generated 
by the electromagnetic field which cuts off the logarithmic interaction HH I22I1 and leads to the absence of any KTB phase 
transition. However, for realistic finite 2D superconducting films where the lateral dimension of the film can be smaller then the 
screening length Rq < X^tf the KTB transition can be restored iSIHt]. This constitutes an intrinsic finite size effect. 

In layered materials, the interlayer coupling modifies the 2D picture and leads to new types of topological defects. If the 
layers are coupled by Josephson coupling (like for many HTSC materials), the vortex-antivortex pairs on the same layer interact 
with each other via a logarithmic term for small distances but they feel a linear confining potential for large distances (see e.g. 
101 and references therein). The vortices in neighboring layers always interact via a linear potential which can couple them by 
forming vortex loops, rings, or vortex "strings" piercing all layers. 

If the layers are coupled by purely magnetic interaction (e.g. in artificially produced superlattices where the Cooper pair 
tunneling between the superconducting layers is suppressed by relatively large insulating layers) the topological defects for a 
system which consists of infinitely many layers are pancake vortices iflollTsIl which undergo a KTB phase transition at Tktb- 
explained e.g. in Ref. jsl], the Josephson coupling can be essentially neglected when the confinement length, i.e. the length scale 
at which the linear confining potential due to the Josephson coupling dominates over the logarithmic interaction due to magnetic 
effects, is pushed beyond the effective screening length for the logarithmic interaction among vortices. This situation is present 
when the tunneling between the superconductinglayers is suppressed by relatively large insulating layers, and a proposal for a 
experimental realization has recently been given jsjl. For a finite number N of magnetically coupled layers, the Pearl type vortex 
stack 12] is broken up into a number of coupled pancake vortices of fractional flux jsl fTsiflsll . and this configuration undergoes a 
KTB-type phase transition at a layer-dependent temperature T^tb ~ ^ktb (1 ~ N~^) which is connected with the dissociation 
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of the stack. This result has been obtained on the basis of the entropy method first introduced in the ground-breaking work 
Recently, a real space renormalization-group (RG) analysis of the case N ^ 2 has been performed in Ref. [5j| using the dilute 
gas approximation. A priori,it appears to be rather difficult to generalize this RG analysis for N > 2 layers. 

In general, the Ginzburg-Landau (GL) theory (23!] provides us with a good theoretical framework in which to investigate the 
vortex dynamics in thin films and in layered materials. Several equivalent models, like field-theoretical, statistical spin models 
and a gas of topological defects have also been used to consider the vortex properties of films and layered systems. The 2D-GL, 
2D-XY and the 2D Coulomb gas models (see e.g. U 0] and references therein) are considered as the appropriate theoretical 
background for the vortex dynamics of superfluid films. The field theoretical counterpart is the 2D sine-Gordon (SG) model [ 2% . 
Both kinds of these models belong to the same universality class and produce the KTB phase transition. For superconducting 
films one has to consider the 2D-GL model in the presence of electromagnetic interactions fi}] or the equivalent gas of topological 
excitations, the 2D Yukawa gas ll2lll . The corresponding field theory is the 2D-SG model with an explicit mass term, the massive 
2D-SG model ll2i|] . 

For Josephson-coupled layered superconductors in the case of very large anisotropy one should investigate the layered GL 
model including the Josephson coupling between the layers 101 (i.e. the Lawrence-Doniach model ll25ll '). In case of not too 
large anisotropy on can use the anisotropic, continuous GL theory which can be mapped onto the isotropic GL model 

by an appropriate rescaling method 112711 . The corresponding spin model is the 3D-XY model luSll and the equivalent gases of 
topological excitations are the layered vortex ll29ll or vortex-loop ll28ll gases. There are attempts in the literature to construct the 
field theoretical countpart of the isotropic model jsot]. In case of strong anisotropy, the layered sine-Gordon (SG) model llsill 
has been proposed as a candidate model where the interlayer interaction between the topological defects has been described by 
a mass matrix which couples the SG fields 

1 J 

n=l 

where ip = {(pi, . . . , (Pn) and (n=l,...,N) are one component scalar fields. Recently, we showed in Ref. Ol] that the layered 
SG model with the above mass matrix is not apropriate for the description of vortex dynamics of Josephson coupled layered 
superconductors . 

In case of purely magnetically coupled layered systems, the layered GL model has to be used but excluding the Josephson 
coupling. Althou gh th e interaction potentials between the topological defects of magnetically coupled layered systems are given 
in Refs. llsilsl fl^llSll . no field theoretical model has been proposed for the description of vortex dynamics in a finite system of 
magnetically coupled superconductors. 

Here, our aim is to open a new platform for considering the vortex dynamics of magnetically coupled layered systems by 
constructing a multi-layer sine-Gordon (MLSG) type field theoretical model where the two-dimensional sine-Gordon (2DSG) 
fields characterizing the layers are coupled by an appropriate general mass matrix, 

\n=l 

By the exact mapping of the MLSG model onto an equivalent gas of topological defects, we recover the interaction potential 
given in Refs. |5, 8, 12, jLSJ and, hence, prove the applicability of the model. We analyse the phase structure of the MLSG model 
by a differential renormalization group (RG) method performed in momentum space, which is in general easier to perform 
than that in real space, and determine the layer-dependence of T^xb- field theoretical RG approach, the RG flow can be 

calculated in one step for an arbitrary number of layers, and the study of the intrinsic finite size effect of thin film superconductors 
10, mi] and of finite layered systems is facilitated. 

This paper is organized as follows. In Sec. [Ill we define the multi-layer sine-Gordon model and show by its exact mapping onto 
the equivalent gas of topological excitations that it is suitable to describe the vortex dominated properties of magnetically coupled 
layered superconductors. In Sec. |llll a renormalization group analysis of the multi-layer sine-Gordon model is performed 
within the framework of the Wegner-Houghton renormalization group method, in momentum space for general A^, and with a 
solution that spans the entire domain from the ultraviolet (UV) to the infrared (IR). The layer-number dependence of the critical 
temperature of the multi-layer sine-Gordon model is determined by using the mass-corrected Unearized RG flow. Conclusions 
are reserved for Sec.lIVI 

II. MULTI-LAYER SINE-GORDON MODEL 

The multi-layer sine-Gordon (MLSG) model consists of A^ coupled two-dimensional sine-Gordon (2D-SG) models of identi- 
cal "frequency" b, each of which corresponds to a single layer described by the scalar fields (pn {n = 1,2,..., A^). Its Euclidean 
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bare action (we imply here the sum over /i = 1, 2) 
contains the interaction terms 



N 



Vn C0s(6(p„) 



(1) 



(2) 



with the 0{N) muhiplet (p = {ipi, . . . , (Pn). We can choose the fugacities y„ > without loss of generality, ensuring that the 
zero-field configuration is a local minimum of the action (see Chap. 31 of Ref. ll32ll ). The mass-matrix describes the interaction 
between the layers and is chosen here to be of the form 



N 



(3) 



where G is the strength of the interlayer interactions, and the a„ are free parameters. As will be explained below, any choice 
with a^j = 1 for all ti = 1, . . . , reproduces exactly the same layer-dependence of r^^B found in Refs. JHHt]. In this case, 
the layers can be assumed to be equivalent and, as a consequence, the fugacity y„ = y for n = 1,2, . . . , N. The most obvious 
choice fulfilling af^ ~ 1, namely a„ = 1 for all n = 1, . . . , iV, reproduces the interlayer interaction between pancake vortices 
given, e.g., in Eq. (89) of Ref. [150, and we will restrict our attention to this choice in the following. 

The MLSG model has a discrete symmetry under the shift of the field variable ip ^ ip + A with A = (Zi27r/fe, . . . , iN^ir/b) 

where the "last" integer 1^ = — '^n=i fixed but all the other integers Z„ (?i = 1, . . . , — 1) can be chosen freely (to see 
this, one just diagonalizes the mass-matrix). The single non-vanishing mass eigenvalue is Mjv = V NG, and hence the model 
possesses A^ — 1 massless 2D-SG fields and a single massive 2D-SG field. After the diagonalization of the mass matrix by a 
suitable rotation of the fields, the model thus is invariant under the independent separate shifts of A^ — 1 massless fields, but the 
explicit mass term of the single massive mode breaks the periodicity in the "massive" direction of the A^-dimensional internal 
space. 

One crucial observation is that the partition function of the MLSG model, whose path-integral formulation reads 



Af / V[p]exp{~Sy]), 



(4) 



can be identically rewritten in terms of an e quiv alent gas of topological excitations (vortices), whose interaction potentials 
are exactly equivalent to those of Refs. iH [sl Il2ll . This finding constitutes a generalization of known connections of the d- 



dimensional globally neutral Coulomb gas and the d-dimensional sine-Gordon model, as discussed in Chap. 32 of Ref. 132 
and can be seen as follows. In Eq. ([T]), one artifically introduces the vectors /„ = . . . , S^n) as projection operators to 



rewrite X]ri=i cos(5i^„) = X]ri=i cos(fe/ ip), one expands the periodic piece of the partition function (IHi in a Taylor series, 
and one introduces the integer-valued charges cTq, = ±1 of the topological defects which are subject to the neutrality condition 
X]q=i '^a = 0. This leads to the intermediate result. 



(Ti ,...,(T(^— ±1 



(5) 



X Jv[ip]cxp - J(fr^(p^ {-d'^ + M^)ip + ibp^ ip 



where 9^ = 9^9^, and 



a=l 



(6) 



We have thus placed the 2;^ vortices, labeled by the index i, onto the A^ layers, with vortex i being placed onto the layer ni. The 
Gaussian integration in Eq. (|5]l can now be performed easily, and the inversion of the matrix —d^ + can be accomplished by 
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going to momentum space. Via a subsequent back-transformation to coordinate space, we finally arrive at the result 

■.2U f 2v N 



(7) 



Vi— 1 rii — \ 



cri....,cri/ — ±1 



exp 



~ ^ 



a, 7 — 1 



where 5nm represents the Kronecker-delta. Equation (O imphes that the parameter 6^ in Eq. (O can naturally be identified as 
being proportional to the inverse of the temperature of the gas, oc T~^. The potentials Aa-y ^ A{fa, t^) and Ba-y = B{ra^r^) 
are the intralayer and interlayer interaction potentials, respectively. They read 



2tt N \ a 



1 1 



1 

" 2tc 



where r^-^ = \ra — and 



' a7 



Aoff / \ Aoff 



[Va-i < Aoff) 
[Va^ > Aeff) 



(8a) 



271- Af \^ Aott 



(r^^ < Aoff) 

{rai > Aoff) 



(8b) 



isTo (r) stands for the modified Bessel function of the second kind, a is the lattice spacing which serves as an UV cutoff and an 
effective screening length Aoff is introduced which is related inversely to the non-zero mass eigenvalue of the mass matrix (O, 
A~ff = Mfq = V NG. The relation Ko{r) = — ln(r) + ln2 — 7e + 0{r) has been used in the derivation of the asymptotic 
short- and long-range forms in Eqs. (ISal i and dSbl i. and only the leading logarithmic terms are indicated (7e = 0.577216 ... is 
Euler's constant). 

The interactionpotentials ([Sj have the same asymptotic behavior as the vortices of magnetic allycoupled superconducting 
layers jsl [III [3] [for the intralayer and interlayer interactions see Eqs. (86) and (89) of Ref. lITsll . under the substitution 
Ko = Kg/N}. This observation shows that the MLSG field theory is suitable to describe the vortex dynamics in magnetically 
coupled layered systems. A few remarks are now in order, (i) The prefactor [N — 1)/N appearing in the intralayer interaction 
indicates the existence of vortices with fractional flux in the MLSG model, (ii) For small distances r <C Aoff, the interlayer 
interaction B disappears and the intralayer potential A has the same logarithmic behaviour with full flux as that of the pure 2D- 
SG model (which belongs to the same universality class as the 2D-XY model and the 2D Coulomb gas). Therefore, the MLSG 
model for small distances behaves as an uncoupled system of 2D-SG models, (iii) For the case = 1, there exists no interlayer 
interaction, and the intralayer potential is logarithmic for small distances and vanishes for large distances. Consequently, there 
are always free, non-interacting vortices in the model which push the KTB transition temperature to zero. The MLSG model 
for a single layer reduces to the massive 2D-SG model discussed in Refs. (0, [iM HH HSl where the periodicity in the internal 
space is broken and the KTB transition is absent, (iv) In the bulk limit N ^ oo, the effective screening length and the interlayer 
interaction disappear (Aoff 0, Bq,^ 0), and the intralayer potential has a logarithmic behaviour with full flux, thus the 

MLSG model predicts the same behaviour as that of the pure 2DSG model with T^^b ~ ^ktb- Alternatively, one may observe 
that for N ^ oo, the effect of the infinitely many zero-mass modes dominates over the effect of the single remaining massive 
mode entirely, leading to a constant limit for the transition temperature as — *■ oo. 

For N = 2 layers, the MLSG model [with the choice a„ = (—1)"^^] has been proposed to describe the vortex properties 
of Josephson coupled layered superconductors Isill . However, the above discussed mapping indicates that any layered sine- 
Gordon model, whatever be the mass matrix, can be mapped onto an equivalent gas of topological excitations, whose interaction 
potentials are determined by the inversion of a two-dimensional propagator of the form —d^+ M^. Any such propagator, upon 
backtransformation to coordinate space, can only lead to a logarithmic behaviour for the vortex interactions at small and large 
distances, and consequently, cannot possibly reproduce the confining linear long-range intralayer interaction given in Eq. (8.42) 
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of Ref. iHl and in Ref. Isill . The candidate fsii] for a mass matrix Lp^ rrf'Lp ~ J "^fjl^i^i — Vi+i)^ has also been discussed 
in Refs. Illl [33l[34ll . This candidate interaction is inspired by a discretization of the anisotropic 3D-SG model Issll . but it cannot 
reproduce the linear confining potential needed for the description of the Josephson-coupled case HI]. The layer-dependent 
transition temperature of this model is Tc oc N^^ and decreases with the number of layers, and for general N, the mass matrix 
TO^ also leads to different short- and long-range intralayer potentials as compared to Eq. ([8]) and cannot be used for the description 
of magnetically coupled TV-layer systems, either Finally, let us note that a suitable model for the Josephson coupled layered 
system could probably be constructed if the interlayer interaction term is represented by a compact variable, i.e., one couple the 
phase (compact) fields between the 2D planes HtIi and not the the dual fields. 



III. RG ANALYSIS OF THE MULTI-LAYER SINE-GORDON MODEL 



The above statements on the MLSG model are based on the bare action where the coupling parameters of the theory are fixed. 
However, only a rigorous RG analysis enables one to construct the phase diagram in a reliable manner. For N ~ 2 layers, the 
phase structure and the vortex properties of the magnetically coupled layered system have already been considered with a real 
space RG approach |5] using a two-stage procedure, and a momentum space RG method jsifl on the basis of the dilute gas 
approximation has also been used. Here, we apply a generalized multi-layer, multi-field Wegner-Houghton (WH) RG analysis 
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FIG. 1 : In the left panels, the mass-corrected scaling [see Eq. (TTJ] of the dimensionless Fourier amplitude 
y of the MLSG model for N = \ (top) and for A'^ = 2 (bottom) layers is represented graphically for \? = 
47r, Stt, 127r, IGtt, 207r (from top to bottom on each panels, see the dashed curves). We use G = 0.0001 in 
order to have the UV and IR regimes conveniently located on the plots, which start at the UV scale A = 1. 
The dotted line is the extrapolation of the UV {k 3> A/jv) scaling to the IR {k <^ Mn) region. For N = 1 
layers, y is always relevant (~ k^^) in the IR. For N = 2 layers, y is relevant for < IGvr in the IR and 
irrelevant for > IGtt. Thus, the 2-layer MLSG model undergoes a KTB type phase transition at b^ — IGvr. 
In general, the KTB transition temperature of the MLSG model is layer-dependent Tj^^b = (1 ^ -'V~^)Tktb- 
If the system has a finite volume (R < oo), the thermodynamic limit cannot be taken automatically and, as a 
simple realization of the finite size effect, a momentum scale fcmin ~ 1 /R appears in the model. For R < X^h 
(i.e. fcmin > V NG = Mn), the phase structure of the MLSG model is determined by the UV scaling which 
predicts a KTB type phase transition at b^ — 8tv for any number of layers. 



developed by us for the layered SG type models to the MLSG model with an arbitrary numbers of layers. 

In the construction of the WH-RG equation, the blocking transformations ifssll are realized by a successive elimination of the 
field fluctuations in the direction of decreasing momenta, in infinitesimal momentum shells, about the moving sharp momentum 
cutoff k (see Ref. fs^ ). The physical effects of the eliminated modes are transferred to the scale-dependences of the coupling 
constants [e.g., y = y{k)]. The WH-RG equation in the local potential approximation (LPA) for the MLSG model with N layers 
reads 



i2 + kdk) Vk = -^In 
An 



dot 



5ij + d^.d^.Vk 



(9) 
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where we have defined the dimensionless blocked potential as 14 = fc ^ T4. We make the following ansatz for the blocked 
potential, 

yk = lGk(^^^,}j +Uki^u---VN), (10) 

where the scale-dependence is encoded in the dimensionless coupling constants y{k) and G{k) which are all related to their 
dimensionful (no tilde) counterparts by a relative factor k~^. Inserting the ansatz (fTOb into Eq. (|9|l, the right hand side becomes 
periodic, while the left-hand side contains both periodic and non-periodic parts Il34i 13611 . 

In order to go beyond the dilute-gas approximation, we calculate a mass-corrected UV approximation of Eq. ^ by expanding 
the logarithm of the determinant in the right hand side of Eq. ^ in powers of the periodic part of the blocked potential. Because 
this procedure has been discussed at length in Refs. Ii33„.34>.36il , we immediately state the result [cf. Eq. (43) of Ref. i36ll 1. 

(JV-l)b^ „ 

yik)^y{A) , - , (11) 



A^ + NGJ \A 

with the initial value y{A) at the UV cutoff fc = A. Let us note that in our RG approach the dimensionful G and are scale- 
independent constants. We can immediately read off from Eq. (fTTI ) the critical value 6^ = 87r/(l — N^^) and the corresponding 

KTB temperature T^^b ^ = ^ktb(1 ^ N^^). The fugacity y is irrelevant (decreasing) for > b^ and relevant 
(increasing) for b^ < b^ for decreasing scale k (see Fig.[T]i. Our RG approach provides a consistent scheme to calculate higher 
order corrections to the linearization in the periodic part of the blocked potential, which is equivalent to higher-order corrections 
to the dilute-gas approximation. For = 1, the mass-corrected UV scaling law ( fTTI ). obtained for the massive SG model, 
recovers the scaling obtained in Refs. ll2lll40ll (no phase transition). 



IV. CONCLUSION AND SUMMARY 



In conclusion, we propose the multi-layer sine-Gordon (MLSG) Lagrangian as a quantum field theoretical model for the 
vortex properties of magnetically coupled layered superconductors. Note that the MLSG model cannot be assumed to belong to 
the same universality class as the layered Ginzburg-Landau model J]]], which entails a discretization of the Ginzburg-Landau 
model in one of the spatial directions. The mapping of the MLSG model onto the gas of topological defects is used to clarify the 
suitability of the MLSG model to magnetically coupled layered systems. We investigate the scaling laws for the MLSG model 
using a functional formulation of the Wegner-Houghton RG approach in the local potential approximation. The linearization of 
the RG flow in the periodic part of the blocked potential (and not in the full potential) enables us to incorporate the effect of the 
interlayer interaction into the mass-corrected UV scaling laws, which improve the dilute gas approximation. The mass-corrected 
Wegner-Houghton UV scaling laws indicate that for general interlayer interactions of the type of Eqs. (O, one finds two phases 
separated by the critical value b^ = Stt/ (1 — N^^), where N is the number of layers. This determines the layer-dependence of 
the KTB transition temperature Tj^^g = ^^ktb ^ N^-^) in full agreement with Refs. OKI}]. Perhaps, further investigations of 
the MLSG model (e.g., beyond the local potential approximation) and other generalizations of the momentum-space RG studies 
presented here could enrich our understanding of the layered structures. 
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